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ABSTRACT 

A combined sample of 79 high and low redshift supernovae la (SNe) is used to set 
constraints on the degree of anisotropy in the Universe out to z ~ 1. First we derive the 
global most probable values of matter density CIm, the cosmological constant and 
the Hubble constant Ho, and find them to be consistent with the published results from 
the two data sets of Riess et al. 1998 (R98) and Perlmutter et al. 1999 (P99). We then 
examine the Hubble diagram (HD, i.e., the luminosity-redshift relation) in different 
directions on the sky by utilising spherical harmonic expansion. In particular, via the 
analysis of the dipole anisotropy, we divide the sky into the two hemispheres that yield 
the most discrepant of the three cosmological parameters, and the scatter Xhd m eacn 
case. The most discrepant values roughly move along the locus —451m + 3J7a = 1 
(cf. P99), but by no more than A w 2.5 along this line. For a perfect FRW universe, 
Monte Carlo realizations that mimic the current set of SNe yield values higher than 
the measured A in ~ 1/5 of the cases. We discuss implications for the validity of the 
Cosmological Principle, and possible calibration problems in the SNe data sets. 

Key words: cosmology: miscellaneous - cosmology: observations - cosmology: theory 
- supernovae:general 



^ ; 1 INTRODUCTION 



1994), the derivation of a n independent rest frame from mul- 



The va idity of the Cosmological Principle and the isotropy 
it impli es gained much credibility in recent years. The small 
fluctuations in the CMB (AT/T ~ 1CT 5 on angular scale ~ 



tiple image lens systems (Kochanek, Kolatt, & Bartelmann 



1996), and Fa raday rotation s ignature due to anisotropic 



10°) provide the strongest evidence that the universe can be 
well approximated by the FRW metric on scales larger than 



magnetic field ( Kronberg 1976 Vallee 1990 
ston 19971). 



Nodland & Ral- 



The recent use of SNe as distance indicators (Phillips 



~ lOOU/i MVlpc (e.g., Peebles 1993; Wu, Lahav, & Rees 
1999) 

On smaller scales (~ 100 ft -1 Mpc) bulk flows of the 
order v/c ~ 10 -3 indicate that this isotropy breaks down. 
This is also manifested by significant correlation functions of 
galaxies and clusters on large scales, and structures like the 
Supergalactic Plane and the Great Attractor. The transition 
scale to isotropy and homogeneity is still poorly known, and 
so is the convergence of the acceleration vector of the Local 
Group with respect to the CMB. It is therefore important 
to quantify the degree of homogeneity and isotropy as func- 
tion of scale. Traditionally this was done by searching for 
anisotropy in the distribution of radio sources and back - 
ground radiations ( Nan fc Cai 199tj ; Evans 1992 ; Webster ). 
Several new methods have been suggested to test isotropy 
and homogeneity on redshift scales of z m 0.1 — 5, such as 



1992 ; Perlmutter et al. 1995 



Riess, Press, & Kirshner 199c) 



opened a new opportunity for accurate measurements of 
anisotropy on cosmological scales that previously have not 
been accessible. So far the SNe have been used in order to 
constrain the Hubble constant Ho from a nearby sample and 
combinations of the matter density Qm and the cosmologi- 
cal constant J1a utilizing SNe at moderate (J; 0.3) and high 
(~ 1) redshifts. In the future, SNe samples over a wider 
redshift range will provide separate estimates for the two 
parameters. It is important to establish the 'universality' of 
the measurements of cosmological parameters from SN, as 
they are commonly used in joint analysis with other probes 
s uch as the CMB , cluster abundance and peculiar velocities 
(|Efstathiou 1999|; [Efstathiou et al. 1999j ; |Bridle et al. 199S| ; 



Bridle et al. 



Tegmark 1999) 



measurements of in situ CMB temperature (Songaila et al 



Assuming a FRW cosmology, a forth measure can be de- 
duced from the 'Hubble diagram' (HD; i.e., the luminosity 
- redshift relation) , the Xhd measure for the best fit model. 
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For a perfect distance indicator this measure indicates devi- 
ations of the local potential (i.e., at the location of the SN) 
from a pure FRW geometry. However, in the real universe 
the deviations can also be due to other sources: 

• Intrinsic (astrophysical) scatter in the SN luminosity- 
light curve relation. 

• Scatter due to the location of the SN within the host 
galaxy & the galaxy type. 

• Scatter due to dust absorption in the host galaxy, in 
the intergalactic medium and in our Galaxy. 

• Gravitational lensing along the l.o.s. to the SN (e.g., 
an overdensity along the l.o.s. will enhance the apparent 
luminosity of a SN). 

Here we explicitly assume that there is no evolution 
with redshift in the luminosity-light curve relation. Fortu- 
nately, most of the abovementioned effects are on the scale of 
the host galaxy, so with large enough sample they would be 
averaged out in the calculation of large scale anisotropies. 
On the other hand, one should worry about 'anisotropies' 
which are simply due to poor matching of different data 
sets that sample different portions of the sky, or large angu- 
lar effects due to Galactic extinction. 

We also note that some of these effects above might be 
correlated with other measurements, e.g. if the scatter Xhd 
detected in SN Hubble diagram is affected by fluctuations 
in the potential, then it would be correlated with Integrated 
SW (or Rees-Schiama) effect in the CMB fluctuations. 

The outline of this paper is as follows, in ^ we present 
the unified data set we will be using for the isotropy analy- 
sis. The results for cosmological parameters from the entire 
sample are presented in the anisotropy measurement is 
discussed in j^, and put in a probabilistic context in We 
conclude our results in 



2 THE UNIFIED DATA SET 

An ideal data set of SNe for the goals we have put forward 
in the introduction would be a whole-sky homogeneous cov- 
erage at various redshifts of SNe. Since such an optimal set 
does not exist, the closest data set would be the amalgama- 
tion of the two existing, published data sets. 

We unify the samples of the S upernova Cosmology 
Project (SCP) ( [Pcrlmutter et al. 1999]) and that of t he High- 
z Supernova search team (HZS) (Riess et al. 1998). These 
include also the data from low redshift of the Calan-Tololo 
survey (Hamuy et al. 1996). The two groups have different 



strategy and different nomenclature for the minimization 
problem by which the cosmological parameters are derived. 
We have brought the SCP data to comply with the language 
of the HZS team 

For each SNe we list its (i) cz in the CMB frame, (ii) 



the distance modulus \i 



j^j fiducial 



(iii) errors for 



these two quantities, (iv) Galactic I and b. For the SCP data 
the fiducial magnitude, Mg duclal (cf. P99), is obtained by 
comparison of the 18 overlapping low redshift SNe as anal- 
ysed by the two groups, and equating the distance modulus 
of R98 (table 10) to ml" 1 of P99 (table 2). This procedure is 
repeated twice, since Riess et al. provide two ways to calcu- 
late the distance moduli, "Mufti Light Curve Shapes" (LCS) 
and "Template" . Errors are taken from the tables and a least 




TtrLCS 
* TEMPLATE 

Figure 1. The sky distribution in Galactic coordinates of the 79 
SNe composing the unified sample. The point size is proportional 
to (1 + of the SNe. Also shown are the (positive) directions 

that maximize the Qm and Oa dipoles using the two methods (cf. 
§^), and the CMB dipole direction in the Local Group rest-frame 
as measured by COBE. 



square minimization is performed in order to obtain the two 



best fit values of y[^ duclal Q f pap, (and to recover the Hub- 
ble constant dependence they omitted in their calculation). 

' + 5 log H = -19.322, -19.453 



fiducial 



The two values are M- 
with x 2 /d-0-f of 0.952 and 0.763 for the LCS method and the 
TEMPLATE method respectively. The value of M f B iducial is 
degenerated with Ho, so different Ho calibrations in the two 
samples get "absorbed" in the value for Mg duclal . The uni- 
fied sample consists of 79 SNe altogether, after the exclusion 
of 6 SNe from P99 (taking their "model C" version) and in- 
cluding the snap-shot survey from R98 along with 1997ck. 
Figure ^ shows the SNe distribution in Galactic coordinates. 

The sky coverage is clearly inhomogeneous: the SNe de- 
ficiency near the Galactic plane is evident and the clustering 
of a few of the observed SNe due to the detection procedure 
is clear. 



3 COSMOLOGICAL PARAMETERS FROM 
THE UNIFIED SAMPLE 

We follow the statistical analysis as described in R98 and 
obtain best values for Ho and probability contours in the 
(Am, Ha) plane after integration (i.e. marginalization) over 
all Ho values and taking into account only physical regions 
in that plane. 

P99 include the error due to redshift measurements and 
peculiar velocities in their magnitude errors, for R98 we 
followed their procedure, set a v = 200 km s _1 for SNe of 
2 < 0.5 and a v = 2500 km s _1 for SNe with z > 0.5, and 
translated to the distance modulus, fi, units according to the 
assumed cosmological model in the likelihood function. Fig- 
ure ^ show the results of the likelihood analysis. The max- 
ima of the likelihood functions are obtained for (Qm,Qa) 
values of (0.40, 0.82) and (0.66, 1.36) for the LCS and TEM- 
PLATE method respectively. The contour lines correspond 
to the 68.3%, 95.4%, and 99.7% confidence levels. 
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4 ANISOTROPY MEASUREMENT 

The natural expansion for anisotropy detection is in spher- 
ical harmonics. The current data are too sparse to allow 
analysis in redshift shells. 

We expand the four two-dimensional parameter 'fields' 
(for Q.m, S1a, Ho, and Xhd) m spherical harmonics. If the 
isotropy assumption is valid we expect deviations from the 
average value to be due to noise, and the angular power 
spectrum should likewise reflect it. This is unless foreground 
effects alter the signal significantly. 

The operational way to calculate the expansion coeffi- 
cients ai m is as follows. 

• Build a random distribution of points ( "mask" ) on the 
sphere. 

• Assign four best-fit parameters to every point based on 
minimization over all SNe within angular radius fmin about 
this grid point. 

• Construct the four residual fields about the global 
mean, i.e., $f — (F — (F))/(F), where F is £Im, ^a, Ho, or 
Xhd- 

• Expand the 5f values as obtained at each grid point in 
Spherical Harmonics up to l ma x = n/jmin, i.e. 

I — Z max m=-\-l 

8 F (e,cj>)= J2 E aTYr - w 

1=0 m=—l 

In order to include more than 2 SNe in each smoothing bin 
(at least two-parameter fit) we obtain 7 ~ 25°, however 
the SNe are not distributed uniformly (cf. Fig. |l|) and thus 
a minimum angular resolution of ~ 60° is imposed. That 
means that for a whole sky coverage the highest significant 
multipole, is / = 3. There are, though regions that are more 
densely covered by SNe data and therefore higher multipoles 
can be assessed as well but at a lower signal-to-noise level. 

In order to account for the Poisson noise contribution 
(and thus to the angular power spectrum in quadrature), we 
run a set of 50 random "masks" and repeat the aim calcu- 



lation each time. For each set of ai m the power spectrum 
coefficients, Ci = (21 + I) -1 _ , |a; m | 2 are computed. 

The angular power spectrum of the Sh is an order of 
magnitude and more smaller than the noise level (Ci/(2n) ~ 
5 x 10~ 5 ), in both methods. Figure ^ shows the angular 
power spectrum, 1(1 + l)Ci / (2-k) , for the other three fields as 
calculated from 50 runs with different random mask points. 
The straight weaker lines show the noise level in each field. 
The Sn fields in both methods show signals that exceed the 
noise level for the dipole (I = 1) and the quadrupole (1 — 2). 
Two factors contribute to the noise level, the discrete num- 
ber of SNe, and the scatter in the luminosity — redshift 
relation. The former is common to both methods (LCS and 
TEMPLATE) and therefore the order of magnitude higher 
noise level for the Sq fields in the LCS method must be 
due to the latter. The TEMPLATE method seems to pro- 
vide smaller errors and a better match between the two data 
sets, as indicated by the lower \ 2 level of the fiducial magni- 
tude calibration (cf. §^|). The S x 2 angular power spectrum 
is similar in shape and magnitude in both methods, and lies 
an order of magnitude to a factor ~ 5 above its noise level. 
This may indicate there exists a true dipole (or quadrupole) 
in this field. From the first multipole of angular power spec- 
trum alone, one cannot deduce what is the dipole direction. 
We therefore turn to look for the direction by other means. 

We search for largest dipole in Qm, ^a, Ho, and Xhd- 
This has been done in two ways : an actual search over the 
sky, dividing the SN population in between two hemispheres, 
and equivalently, by solving a maximization problem of the 
dipole term with respect to (8, (j>) using the computed aim 
coefficients. Both methods yield similar results. We then cal- 
culate the confidence regions for each hemisphere separately, 
and look for statistical consistency (overlapping contours). 
Each test can be applied to each one of the four parameters. 

Figure ^ verify the fact that the current SNe data best 
constrain a linear combination of the cosmological parame- 
ters Om, ^a- In all four panels the likelihood maxima move 
along the line (P99) — 4Qm + 3Qa ~ 1, sometimes with 
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Figure 3. Angular power spectrum coefficients for the two dimensional fields: Sq , Sq, , and S 2 along with the noise level of each 
field (straight weak lines). Shown are results from the LCS method (right) and the TEMPLATE method (left). 



a large distance A between the two maxima for the two 
disjoint hemispheres (quoted on the plots). In three cases 
the contour levels overlap significantly (see next section for 
quantitative evaluation). In the case of the £7a dipole, us- 
ing the LCS method, there is no overlap between the 99.7% 
confidence levels of the two hemispheres. The discrepancy 
stems from the very assymetric distribution of SNe between 
the two hemispheres (59 on one versus 20 on the other) and 
only one SN (1995at) with z > 0.5 in the 20 SNe sample. 
A small error in the distance measurement of this SN, or 
a systematic deviation of it from the average LCS relation 
may cause such a discrepancy as we demonstrate in the next 
section. Elimination of this SN yields a dipole which points 
~ 20° away from the original direction, reduced A value 
of 2.39, and almost full inclusion of the 99.7% confidence 
contour for the larger sample (56 SNe) within the 95.4% 
confidence level of the remaining 22 SNe. 

Note that a different "mixture" of redshift distribution 
to different directions may cause some directions to become 
more sensitive to one parameter. E.g, SNe at z ~ 0.4 — 0.5 
are mostly sensitive to the Q m — combination, as opposed 
to higher weight on Hm as redshift increases (|(9A/<9f2M| > 
\0A/dflA\)- Figure [i] includes the dipole directions (positive) 
of Q in both methods. We observe no coincidence with any 
Galactic or CMB direction, moreover not all dipoles point 
to the same direction. 

The dipole of the 5^2 field points in both methods 

"hd 

toward (I = 80°, b = -20°) with (x H d> = 1.00(1.01) and 
largest difference of 0.60 (0.63) for the LCS (TEMPLATE) 
method. This dipole direction is suspiciously close to the 
Galactic plane. 

One worry is that the detected signal is due to the 
(mis)match between the two data sets. We therefore re- 
peated the computation for each data set separately and 
verified that though the noise level increases, the results as 
drawn from each one of the data sets are consistent with the 
results from the unified set both in magnitude and direction. 



5 DEGREE OF ANISOTROPY 

The results of the last section, regarding the spherical har- 
monic expansion and the various dipole magnitudes, should 
now be put in an expected distribution in order to draw 
conclusions about the degree of anisotropy. 

The hypothesis we are trying to address is that the 
SN data do not falsify the FRW geometry as a reliable de- 
scription of the z ~ 1 Universe. This strategy is more effi- 
cient ,j4ianja£kln3s^ cosmological mod- 



els (Sclerier 2000a; |Celerier 2000b). We therefore compute 



the probability distribution of the dipole magnitudes within 
a FRW universe and confront it with the values obtained for 
the real Universe. 

A simple two dimensional Kolmogorov-Smirnov test to 
falsify the hypothesis that the two contour maps come from 
the same underlying distribution of cosmological parame- 
ters is inadequate here. Since we have used the maximum 
discrepant values in order to obtain the dipole, the two sub- 
samples are not randomly selected and therefore can not be 
confronted in a KS test. 

The probability distribution depends on the actual cos- 
mological values and to a lesser extent on the power spec- 
trum (via the scatter due to potential fluctuations). For a 
self consistency check, the underlying cosmology is taken to 
be the "best fit" cosmological model (§|3|), which we then 
sample by Monte-Carlo simulations. 

To mimic accurately the SN sample, we use the same 
angular locations and redshift values as of the observed sam- 
ple. Luminosity distances, magnitude scatter and peculiar 
velocities are drawn from Gaussian distributions with the 
appropriate observed standard deviation. 

The dipole analysis is repeated for 200 mock catalogs 
of the SN and the maximal dipole magnitude is calculated 
to obtain its distribution for the current sampled SNe. 

Table 1 shows the rejection levels of the hypothesis that 
the Universe up to z ~ 1 can be described by a FRW metric. 
E.g., using the LCS method and the current sample of SNIa 
we expect in 19% of all cases to detect a higher A value for 
Q.M dipole, than the observed one. 
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Figure 4. Confidence regions (same as Fig. Q drawn from 79 SNe of the unified sample in two hemispheres that maximize the fijj 
dipole (up) and £2 a dipole (bottom) using the LCS method (left) and the TEMPLATE method (right). Marked are the number of SNe 
in each hemisphere, the best-fit Ho, and the distance between peak probabilities (A). 



Table 1 



Isotropy rejection levels using A 


Cosmological 




Method 


parameter 


LCS 


TEMPLATE 


H 


33% 


70% 




81% 


79% 




88% 


64% 



6 DISCUSSION 

By the exploitation of the current available SNe data we 
have put constraints on the rejection level of the cosmolog- 
ical principle validity up to z ~ 1. A FRW metric is found 
to be an adequate description of the Universe. In ~ 20% 
of all realizations of such universes, the dipole signature for 



anisotropy in the cosmological parameters Ho, SIm and Q,a 
exceeds the observed one. 

Even though such dipole magnitudes are reasonable in 
the framework of the FRW model, they may be indicative of 
non-cosmological contributions to the angular power spec- 
trum. In §[j] we listed possible such contributions. If indeed 
the Universe up to z ~ 1 is well represented by a FRW 
metric then we can exclude large coherent structures at 
z ^ 0.3. Such are the structures that may lead to dipole and 
quadrupole signatures due to coherent gravitational lensing 
magnification/de-magnification and therefore the latter can 
be excluded as anisotropy contributors. 

That leaves small scale (Galactic) foreground effects to 
be the most likely power contributors. The Galactic disk 
geometry makes the quadrupole the most significant multi- 
pole to be considered, though the solar system offset from 
the Galactic center may bring about a dipole contribution 
as well. In the current sample the quadrupole term is only 
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slightly larger than the noise level and does not allow any 
conclusive results. None of the dipole directions for f2 coin- 
cides with the Galactic plane and thus they are probably not 
correlated with it. Multipoles due to dust extinction may be 
affirmed by multiple expansion of the residual colors after 
extinction correction (i.e., R98 and P99 appendices). 

The one case where two significantly non-overlapping 
confidence regions are found for two hemispheres that max- 
imize the Oa dipole (LCS), is probably due to a single SN 
(1995at) for which the individual errors have been underes- 
timated. This case is an exception since the overall \ 2 val- 
ues for the HD fits are statistically acceptable. Nevertheless, 
this case demonstrates the hazard in the draw of conclusions 
based on a handful of SNe, for which the error in the error 
estimate is uncertain. 

In general, the TEMPLATE method provides a better 
statistical agreement of the data with an FRW model and 
the current SNe data. This is seen from the magnitude match 
(cf. §^), tighter constraints from the combined set, smaller 
noise levels for all multipoles, and smaller A values for the 
Q dipoles. 

We conclude that an isotropic z ~ 1 universe cannot be 
rejected by more than a lcr level based on the current SNe 
data. 
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